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Aaronson and Arkhipov showed that predicting or reproducing the measurement statistics of a 
general linear optics circuit with a single Fock-state input is a classically hard problem. Here we show 
that this problem, known as boson sampling, is as hard as simulating the short time evolution of a 
large but simple spin model with long-range XY interactions. The conditions for this equivalence 
are the same for efficient boson sampling, namely having small number of photons (excitations) 
as compared to the number of modes (spins). This mapping allows efficient implementations of 
boson sampling in small quantum computers and simulators and sheds light on the complexity of 
time-evolution with critical spin models. 


Boson sampling requires (i) an optical circuit with M 
modes, randomly sampled from the Haar measure; (ii) an 
input state with N M photons, with at most one pho¬ 
ton per mode; (iii) photon counters at the output ports 
that post-select events with at most one photon per port. 
Under these conditions, the probability distribution for 
any configuration n S p(ni, 712 ... tim) = iTnP, is 
proportional to the permanent of a complex matrix whose 
computation is #-P hard. This result, combined with 
some reasonable conjectures [T], implies that linear op¬ 
tics and interferometers have computing power that ex¬ 
ceeds that of classical computation, and that the classical 
simulation of random optical circuits with non-classical 
inputs likely involves itself an exponential overhead of 
resources. More recently, boson sampling has been gen¬ 
eralized to consider other input states w , extensions to 
Fourier sampling [3] or trapped ion implementations 
Boson sampling has also been related to practical prob¬ 
lems, such as the prediction of molecular spectra and 
quantum metrology [7]. Finally, there are other quantum 
models, such as circuits of commuting quantum gates [5] 
which also establish potential limits of what can be clas¬ 
sically simulated. 

In this article we prove that boson sampling is equiv¬ 
alent to a many-body problem with spins that interact 
through a long-range, XY coupling and evolve for a very 
short time, of the order of a single hopping or spin-swap 
event. The model involves a bipartite set of input and 
output spins 

M 

+ H.C., (1) 

joined by the (unitary) matrix R. We show that the time 
evolution of an initial state that has only N excited in¬ 
put spins approximates the wavefunction of the boson 
sampling problem after a finite time t = tt 12. All errors 
in this mapping can be assimilated to bunching of exci¬ 
tations in the optical circuit, and the mapping succeeds 
whenever boson sampling actually does. 

Where does Hamiltonian ([^ come from? It arises from 
a setup with two-level systems as photo-emitters and de¬ 
tectors at the entrance and exit of an M-port interferom- 
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FIG. 1. (a) A setup that consists of beam splitters and free 
propagation implements boson sampling if the input state has 
a fixed number of single photons on each port (red wiggles), 
(b) We can regard those photons as arising from the sponta¬ 
neous emission of two-level systems onto the circuit (up spins), 
which after propagation map onto other two-level systems at 
the end, via the unitary matrix R. 


eter [cf. Fig. Efi- An interpretation of this setup in terms 
of earlier works with qubits in photonic waveguides mm 
shows it contains a coherent, photon-mediated interac¬ 
tions of the form Q, accompanied by collective dissipa¬ 
tion terms. In addition of this physical motivation [12) . 
our work introduces a new physical problem that has the 
same complexity as boson sampling, but which can be 
reproduced on state-of-the art quantum simulators and 
small quantum computers. Unlike boson sampling, spin 
sampling may be implemented using error correction, an 
idea that pushes the boundaries of what can be exper¬ 
imentally implemented beyond the limitations of linear 
optical circuits. 


RESULTS 

The boson sampling Hamiltonian 

Before studying the spin model 0, let us develop a 
Hamiltonian for boson sampling. Let R a the unitary 
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transformation implemented by the circuit in Fig. [^. R 
is sampled from U(M) according to the Haar measure, 
and has associated an equivalent Hamiltonian 

M M 

Hbs = ^ +H.C.)+a]aj). (2) 

i,i=i i=i 

for the input and output modes, a and b, of the problem. 
Evolution with this Hamiltonian transforms the initial 
state of boson sampling 

1^(0)) =al---a]v|vac), (3) 

into the N boson sampling superposition 

N 

|(/)(7r/2)) = {-i)^ P |vac), (4) 

i=l j 

with a photon distribution given by the permanents 
|7nP = I (vac|6|"' • • • 6]^^" |</>(7r/2)) p, n* S {0,1}. 

Dilute limit 

The hnal state in Eq. Q contains a non-zero proba¬ 
bility of two or more bosons accumulating in the same 
mode. It can be split as in Fig. 

|(/i) =Q\cj)) + |e), (5) 

where Q \ 4>) is the projection onto states with zero or one 
boson per site and je) contains bunched states. For the 
errors |e) to be eliminated in postselection while main¬ 
taining the efficiency of the sampling, the number of 
modes must be larger than the number of excitations. 
Formally the limit seems to be M ~ log^ (N) , while 
M ~ N'^ is the suspected ratio [1] at which sampling 
becomes efficient, with bounds being tested theoretically 

and experimentally 1111131 - 

Spin sampling 

The assumption of ’diluteness’ of excitations, which is 
needed for the efficient sampling of bosons, not only en¬ 
sures that we have a small probability of boson bunch 
at the end of the beam-splitter dynamics, but also at all 
times. More precisely, the probability of bunching of par¬ 
ticles in the state estimated by ||Q(/'(t)|| 2 i roughly 

grows in time and is bounded by the final postselection 
success probability [cf. Ref. [IS]]. In other words, bo¬ 
son sampling dynamics is efficient only when it samples 
states with at most one boson per mode, the so called 
hard-core-boson (HCB) subspace or the spin space. In 
this situation one would expect that the models 0 and 
0 become equivalent, with the soft-boson corrections 
becoming negligible. 



Hard core bosons 

|(^(0)) = a\ ■ • -ajvlvac) 

|!/'(0)) =(T+---(7+|0i---02m) 

FIG. 2. The distance between the full bosonic state, \4>{t)} 
and the state approximated with spins, \ip{t)), is covered by 
two error vectors: one |5) that lives in the hard-core boson 
space (red dashed), and another one that covers the distance 
between the projected state Q \4>{t)} and the full boson state, 
\(j)} (black dashed). 

Continuing with this line of thought, we now study 
how well the dynamics of the full bosonic system can be 
approximated by the hard-core boson model Q. We re¬ 
gard the spin Hamiltonian as the projection of the full bo¬ 
son sampling model onto the hard-core-boson subspace, 
H = QHbsQ- Using this idea, we show that the boson¬ 
sampling dynamics is reproduced by this spin model at 
short times, with an error that grows with excitation den¬ 
sity, and which is bounded by ||£p. 

Let us assume that \'ip) is a hard-core boson state that 
initially coincides with the starting distribution of the 
boson sampling problem, |?/'(0)) = |(/'(0)). This state 
evolves with the spin model as idt \ip) = QHbsQ W ■ Let 
us introduce the actual sampling error which we make by 
working with spins 

|<5) =Q\(j)) - \tjj). (6) 

A formal solution for this error 

\d{t)) = -t f |e(r)) dr. (7) 

Jo 

Bounds for the different terms in this integral provide us 
with the core result (see supplementary material [15]): 

\\S{t)h<txo(^^y ( 8 ) 

which states that the error probability at t = 7r/2 is negli¬ 
gible when ~ C>(M^U), jn that material we also show 
numerical evidence m that this bound can be at least 
improved to iV ~ Moreover, we show that the 

same bounds apply for the variation distance between 
the probability distributions associated to the quantum 
states Q \ (j)) and \ip), the measure used by Aaronson and 
Arkhipov in Ref. [T|. 

Equation Q shows that the errors due to using a spin 
model for boson sampling feed from bunching events in 
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the original problem, |(5) oc |e), which are prevented by 
the HCB condition. For times short enough, these er¬ 
rors amount to excitations being “back-scattered” to the 
“in” spins. This means that sampling errors can be effi¬ 
ciently postselected in any given realization of these ex¬ 
periments, rejecting measurement outcomes where there 
contain less than N excitations in the ^ spins. In 
this case, what we characterized as an error becomes a 
postselection success probability, Pok = 1 — ||<J|P- 


Quantum computing 

Our previous results suggest an efficient implementa¬ 
tion of boson sampling in a general purpose quantum 
computer using the following algorithm: (i) Prepare a 
quantum register with M + M qubits, encoding the in¬ 
put and output spins, respectively, (ii) Initialize the reg¬ 
ister to the state |^/’(0)) = |Ii,..., I^v, Ojv+i,..., 02 m)- 
(iii) Implement the unitary U = exp(—ii77r/2), where 
H is given by Q. (iv) Measure the quantum register. 
Postselect experiments where the N qubits are at zero, 
recording the resulting state of the output qubits to esti¬ 
mate the sampling probability. Note that step (iii) may 
overcome the accuracy limitations of optical devices with 
the use of error correction, allowing scaling to arbitrar¬ 
ily large problems. It is worth mentioning that while a 
general purpose quantum computer might implement bo¬ 
son sampling via the Schwinger representation of bosons, 
this would require larger number of qubit resources, and 
a greater complexity in implementing beam-splitting op¬ 
erations, while our spin sampling problem requires a 
smaller Hilbert space, and has a natural implementation 
on a small quantum computer. 


Quantum simulation 

We can use a quantum simulator with spins to imple¬ 
ment spin sampling. As a concrete application, let us 
assume that we have a quantum simulator that imple¬ 
ments the Ising model with arbitrary connectivity and 
coupling to a transverse magnetic field 

2M 2M 

^ (9) 

a,b—l a—1 

In the limit of very large transverse magnetic field, 
\B\ ^ II J||, we can map this problem via a rotating wave 
approximation to the Hamiltonian 0 where the coupling 
matrix is Jij+M = J*+m i = (bi = 1, ■ • ■, M). 

The Ising interaction 0 is already present in trapped 
ions quantum simulators with phonon-mediated interac¬ 
tions [16], a setup which has been repeatedly demon¬ 
strated in experiments HlHIS] , even for frustrated models 
[101 HD, extremely large 2D crystals [Hj, and in partic¬ 
ular in the XY limit [13]. Another suitable platform for 


this kind of simulations is the D-Wave machine or equiva¬ 
lent superconducting processors with long-range tunable 
interactions [HHH]. These devices can now randomly 
sample J from a set of unitaries over a graph that is 
a subset of the available connectivity graph. Since the 
number of spins is very large, with over 900 good-quality 
qubits available, we expect that those simulations would 
surpass the complexity of the sampling problems that can 
be modeled in state-of-the art linear optics circuits. 


Complexity theory 

Our mapping of boson sampling to spin evolution 
shows that classically simulating the dynamics of long- 
range interacting spin models at short times has exactly 
the same complexity, which if the conjectures in Ref. [T] 
hold, is #-P hard. More precisely, if spin sampling were 
to be solvable in a classical computer, then we would 
approximate the boson sampling solution with precision 
poly(A^^/'\/M), which would imply a collapse of the poly¬ 
nomial hierarchy [1]. We have thus established a new 
family of problems that are efficiently simulatable in a 
quantum computer but not on a classical one. 

This idea connects to earlier results that relate the dif¬ 
ficulty of classically simulating time-evolution due to very 
fast entanglement growth [271128]. It also does not con¬ 
tradict the fact that free fermionic problems can be ef¬ 
ficiently sampled because model Q only maps to free 
fermions for a subclass of matrices, R, which are tridiag¬ 
onal. 

There are other remarks to be done about our work 
and its place in the existing literature. First of all, it can 
be argued that spins or qubits are the underlying compo¬ 
nents of a quantum computer whose computation will in 
general amount to evolution with an effective Hamilto¬ 
nian. This argument is bogus in that the resulting Hamil¬ 
tonian will, in general, not be physically implementable, 
involving interactions to an arbitrary number of spins and 
distance. Moreover, even if certain models such[2are uni¬ 
versal and may encode quantum computations [1I1I3D], 
the timescales of our result amount to a single hopping 
event, which is scarcely the time to implement a single 
quantum gate and not an arguably complex computation. 
Finally, while at least one work has established connec¬ 
tions between the collapse of the polynomial hierarchy 
and spin models [8], that works builds on the conjecture 
that the complex partition function of a spin model is 
already in #-P, and thus time-evolution of those spin 
models is hard to be approximated, which is instead the 
conclusion of this work. 

Summing up, we have established that boson sampling 
can also be efficiently implemented using spins or qubits 
interacting through a rather straightforward XY Hamil¬ 
tonian. This map opens the door to simulating this prob¬ 
lem with quantum simulators of spin model, of which we 
have offered two examples: trapped ions and supercon¬ 
ducting circuits. Moreover, the same map states that 
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boson sampling can be efficiently simulated in a general 
purpose quantum computer. 

METHODS 
Linear model 

We build new orthogonal modes c] = Rjia\, that 
satisfy the appropriate bosonic commutation relations, 
[Cm,4] = R)m,n = 5n,m- These 

canonical modes almost diagonalize the previous Hamil¬ 
tonian, which becomes a sum of beam-splitter models 
Hbs = The dynamics of this Hamil¬ 

tonian involves a swap of excitations from the normal 
modes Ci into the output modes bi, so that after a time 
t = 7r/2 the initial state ([^ is transformed into 0. 

It is very important to remark that the state \(nt)) can 
at all times be written as 

N 

l'/'(^)) = cos(t)” sin(t)^""' , (10) 

71 — 0 

where \^n,M) i® output state of a Boson-Sampling 
problem with n input excitations in M modes. Thus, 
has the bunching statistics of efficiently sampled 
Boson-Sampling problems, and it is ruled by the boson 
birthday paradox [T]. This in in sharp contrast with 
actual intermediate states in a random interferometer, 
which may contain a lot of bunching before the bosons 
exit the circuit, and it is due to the fact that the model 
that we use to recreate the BS output states, H, does not 
describe those intermediate stages, where the dynamics 
of individual beam-splitters matters. 

A. Spin model bounds 

Equation 0 arises from the simple Schrodinger equa¬ 
tion 

idt\S) = QHbsQ\S) + QHbs Is) ■ (H) 


As explained above, it shows that the errors in approx¬ 
imating the boson sampling with spins result from the 
accumulation of processes that, through a single applica¬ 
tion of HbS) undo a pair of bosons from £, taking this 
vector into the hard-core boson sector. 

We now bound the maximum error probability as an 
integral of two norms. For that we realize that out of 
varepsilon^ QHbs cancels all terms that have more than 
one mode with double occupation. Thus, 

= ||< 5||2 < [ \\QHBsPlbpatrh\\Plbpatr\e{T))\\2dT, 

Jo 

( 12 ) 

where Q is a projector onto HCB states with N particles 
and Pibpair is a projector onto the states with TV — 2 iso¬ 
lated bosons and 1 pair of b bosons on the same site. As 
explained in the supplementary material m, the value 
\\Pibpair \£{t))\\'^ = \\Pibpair \4'M)\\‘^ is the probability 
of finding a single bunched pair in the full bosonic state. 
Combining a similar bound by Arkhipov [T] with the ac¬ 
tual structure of the evolved state, we find 

\\Plbpair |e(T)) II 2 ^ O ^ 

which works provided that N = We have also 

shown m that the operator norm \\QHBsPibpair \\2 is 
strictly smaller than the maximum kinetic energy of N 
bosons in the original model, Hbs, so that 

\\QHBsPlbpa^r\\2 < N. (14) 

Combining both bounds we finally end up with ([^. 

Note that because we only work with our bound for 
times ||H|| 2 t < 7r/2, the Boson-Sampling error je) con¬ 
tains only some multiply occupied sites at the out modes, 
bj. Moroever, because of the construct in 6, the domi¬ 
nant contribution from e to the error consists on emp¬ 
tying a single bunched site bj and using it to refill an 
input boson, ak , leading to a spin configuration where an 
excitation has been “back-scattered” to (Jin,k- 
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Supplementary material 


I. BOSON SAMPLING DYNAMICS 

Let us begin with the beam-splitter Hamiltonian 
Hbs = 

3 

defined in terms of the transformed modes 

M 

c]='^Rjial. (16) 

We need to study how the evolution of an initial state 
with N M excitations 


N 


iV’)=n 4 io). 

(17) 

k=l 


For that we write down the Heisenberg equations 
erators evolving as 0{t) = 

for op- 

±b] = -i[H,b]] = -zcl, 

(18) 


(19) 

which has as solutions 


6](t) = cos(t)6](0) - isin(t)c](0). 

(20) 

c](t) = cos(t)c](0) - jsin(t)6](0). 

(21) 

Inverting the relation (|16[), we recover 


4w = 5I^4-c]w 

(22) 


= cos{t)R*^Rjial{0) - ism{t)R*f,bl{0) 

= cos(t)aj,(0) - isin(t)i?*j,6](0), 

where we implicitly assume summation over repeated in¬ 
dices. Dynamics under Hamiltonian ( [T^ is coherently 
transferring population from the a to the b modes, as in 

N 

W)) = n “ * sin(t)i?*fc6] j |0). (23) 

k^l 

At time t = tt 12 all population is transferred 

N 

l0(7r/2)) = n 4(^/2) |0) : with (24) 

k^l 

M 

4(^/2) = 

i=i 

But in-between, \(l>{t)) may be regarded as a coherent su¬ 
perposition of different boson-sampling instances, ^BS,n 


where only n = 0,1..., bosons participate and fed into 
the M output modes, while N, N — 1... ,0 remain in the 
input modes. In other words, we have 

I'/'W) = ( ) cos(t)^"" sin(t)” , (25) 

n=0 '^^2 

with normalized states 

It is important now to discuss each of these states, 
which we may write as 

l^BS,n) OC 51 4iW2)---4„(V2) X (26) 

{fed} 

x44o)---4«-4o)Io)’ 

and which consist on N — n excitations that stay in the 
input modes (aj(0)) and n excitations that have been 
fully transferred to afc(7r/2), which are linear combina¬ 
tions of the 6fc(0) modes. In other words, each of these 
instances ICss.n) represent themselves the outcome of a 
Boson Sampling experiments with n input excitations 
and they will have the properties and statistics of bunch¬ 
ing events of any boson sampling problem of such size [T] , 
a fact that we will use for proving bounds on the distance 
between \(j>{t)) and the hard-core-boson subspace at all 
times. 


A. Comparison with Optics 

It is very important to understand while the bunch¬ 
ing statistics of |0(t)) does not contradict the behavior 
of actual optical circuits. In this work we are only study¬ 
ing the unitary transformation implemented by a boson 
sampling circuit. That transformation is generated by 
our toy Hamiltonian, Hbs 

Wbs ■= eM-^HBs7r/2), (27) 

but intermediate stages have no other physical reality 
than being an aid in proving our other results regard¬ 
ing the separation between boson and spin transforma¬ 
tions. In particular, the optical transformation Wbs will 
in general be implemented as a sequence of elementary 
transformations —beam splitters and phase shifters—, 

Wbs = \{W,, (28) 

i 

and intermediate stages of those transformations will, 
in absolute generality, lead to highly bunch and also 
highly entangled states which may be useful for other 
purposes. However, the final state Wbs |?^>(0)) can not be 
highly bunched, or otherwise it would not fit the frame¬ 
work of boson-sampling (And indeed it is not, following 
Arkhipov’s boson birthday paradox). 
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II. BUNCHING BOUNDS 


If we want Boson Sampling to be efficient, we need to 
impose that the number of bunching events in \4>{'k/2)) 
remains small with increasing problem size. Such prop¬ 
erty is guaranteed on average by the random unitaries 
U sampled with the Haar measure, as explained by 
Arkhipov and Kuperberg in the boson-birthday paradox 
paper [1]. Below we will use the fact that the number 
of bunching events in |0(7r/2)) = \^bs,n) indeed upper- 
bounds the number of bunches in each of its constituents, 
\^BS,n<N), and use this idea to draw conclusions on the 
distance between the true Boson Sampling problem and 
the HCB spin model. 

We cannot sufficiently stress the fact that the number 
of bunching events in |0(t)) is not related to the num¬ 
ber of bunching events in the intermediate stages of a 
linear optics circuit. In order to implement a boson sam¬ 
pler one has to combine beam splitters that at differ¬ 
ent stages of the circuit cause the accumulation of the 
bosons. However, while those intermediate states in the 
construct are essential to reach the final Boson Sampling 
states, none of those intermediate states belongs 

to the family of states at the output of the circuit, , 

which must have a low bunching probability (and which 
do, as shown by Ref. m)- 


A. One-bunching events 


For the purposes of bounding the error from the spin¬ 
sampling model, we need to bound the part of the error 
je) that contains a single pair of bosons on the same site, 
on top of a background of singly occupied and empty 
states. We have labeled that component HPipair^lli- 
However, as discussed in Ref. |T], bounding that probabil¬ 
ity is harder than bounding the probability phcb{N, M) 
of having no bunching event in a state with N bosons in 
M modes, distributed according to the random matrices 
Rji. This probability is 


Piicb{N,M) 


N 


n 

a=0 


M -a 


M + a 




(29) 


for dilute systems N = We now use (i) that 

the state \4>{t)) in Eq. ( p^ is made of a superposition of 
states with n = 0,1,... N bosons distributed through the 
M modes, (ii) that due to the randomness of R, each of 
these components shares the same statistical properties 
of the boson-sampling states [T] , (iii) the probability dis¬ 
tribution pbcb{N, M) is monotonously decreasing with 


N. Using Eq. (25) and this idea we arrive at 

\\QHt)\\l - X! ( ) cos(t)^(^"") sin(^)2>HCB(^^, M) 

n=0 '^^7 


> ^ cos(t)2(^-")sin(t)2>HCB(Af,M) 

n=0 

= (cos(t)^ +sm{ty)^p}iCB{N,M) 

= PncBiN,M). (30) 


Using the fact that Q \(j)) and |e) are orthogonal and thus 
Il<^ll 2 = llO'/'lli + Ikilij can find a very loose bound for 
the error probability of single bunching events 

\\Plbpa^re\\l < \\eit)\\l <1 -pbcb{N,M). (31) 


Note that this bound can be translated into an upper 
bound of 0{N'^/M) using the fact that the exponential 
falls faster than 1 — N‘^/M. 


III. HCB OPERATOR BOUND 

In addition to bounding the error vector, we also need 
to bound the norm of an operator that brings back pop¬ 
ulation from the error subspace into the hard-core-boson 
subspace. Because \\Pibpair £\\2 is already rather small, we 
can afford a loose bound for the operator \\QHBsPibpair\\, 
which is the other part of the integral. The argument is 
basically as follows. Eirst, we notice that all operators 
in the product, Q,Hbs and Pibpain commute with the 
total number of particles, which in our problem is exactly 
N. We can thus study the restrictions of these operators 
to this sector, which we denote as PnOPn for each op¬ 
erator, where P/v is the projector onto the space with N 
particles. We then realize that ||AH ||2 < ||A|| 2 ||i ?||2 and 
since the projectors have norm 1, 

WQHBsPlbpairh = || || 2 (32) 

^ llQIbll^Af^^BS^Wlbll-Plbpairlb 

= \\PnHbsPn\\2 

Notice now that PnHbsPn is just the Hamiltonian of N 
free bosons, without hard-core restrictions of any kind. 
In other words, it is the restriction of 

HBS = Y.(blck+li.c.) (33) 

k 

to a situation where cj(,Cfc -|- b\hk = N. We intro¬ 
duce superposition modes, ak± = {ck ± bk)/V^, and 
diagonalize 

Hbs = ^(al+afe+ - al_ak-), (34) 

k 

where the constraint is the same Oi\j^ak+Poi\-^k- = 
N. Since the largest eigenvalues (in modulus) are ob¬ 
tained by filling N of these normal modes with the same 
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FIG. 3. Numerical estimates of the norm \\QHBsPibpair \\2 as 
a function of the number of bosonic modes, M, for different 
number of excitations, N. 


frequency sign, we have 

\\QHBsPlbpa^r\\2 < \\PnHbsPn\\2 = N. (35) 

Note that this proof does not make use of any prop¬ 
erties of H such as the fact that it is built from random 
matrices. As explained in the body of the letter, if we 
sample QHBsPibpair randomly with the Haar measure 
and average the resulting norms, the bound seems closer 
to 0{y/N). 

We have strong evidence that this bound can be sig- 
nihcantly improved using the properties of random ma¬ 
trices Rij and the structure of QHbsS- In particu¬ 
lar, we have numerical evidence that the average norm 
over the Haar measure is \\QHBsPibpair \\2 oc. 
which improves the requirement for efficient spin sam¬ 
pling N ~ Fig. shows the average and stan¬ 

dard deviation of the operator norm obtained by sam¬ 
pling random bosonic circuits with N — 2 — % particles 
in M = 7 — 60 modes, creating random unitaries ac¬ 
cording to the Haar measure and estimating the norm of 
the operator QHBsPupair with a sparse singular value 
solver. Note how, despite the moderate sample size (200 
random matrices for each size) the standard deviation is 
extremely small, indicating the low probability of large 
errors and the efficiency of the sampling. 



IV. VARIATION DISTANCE 


the variation distance between probability distributions 
(this is, 1-norm). However, our proof above can be writ¬ 
ten in a similar was as the one given by Aaronson and 
Arkhipov, that is 

\pi - P 2 \i - P 2 (n)I, (36) 

n 

which represents total difference between probabilities for 
all configurations n = (ni,..., um) of the occupations at 
the output ports. In our model, the probability distribu¬ 
tion associated to boson sampling would be 

Pi(n) = |(n|0(t))|^ =: |(/)(n)|2, (37) 

where if we focus on events with Ui € {0,1}, we can 
replace (j) with Qcj). The corresponding probability for 
the spin model would be 

P2(n) = |(n|V'(t))|^ =: |V’(n)P, (38) 


Using the above expressions for the probability distri¬ 
butions, we can write down the following identities for 
the total variation distance ( [M| ). 

^ |pi(n)-p 2 (n)| = ^||V'(n)|^- |(/>(n)H 

n n 

n 

= ^ X + ('^(n)IV’(n) + <('(n))]|, 


where |5(n)) = |'!/'(n) — ^!)(n)). Hence, it follows that 
IIV'(n)P - |<(>(n)n = X |IIe((V'(n) + ^(n)|5(n)))| 

n 


= '^\^e{2(l)*{n)S{n) + 6*{n)S{n)\ 

n 

< 2Xl<^(n)||<5(n)|+Xl^(n)P (40) 

n n 

= 2 ||<^|| 2||^||2 + || 5 |!i 


Since the boson sampling wavefunction is normalized, 
||(()||2 = 1, and 11^112 < 1 we finally get the next tight 
bound for the variation distance 


Throughout this manuscript, we have found bounds 
according to the 2-norm, in contrast to Aaronson and 
Arkhipov work, whose results are expressed in terms of 


\Pi 


P2\<m\2 = o 



(41) 
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